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Abstract 

In this paper we study the properties of vortices in two dimensional quantum 
antiferromagnets with spin magmitude on a square lattice within the frame- 
work of Schwinger-boson mean-field theory. Based on a continuum descrip- 
tion, we show that vortices are stable topological excitations in the disordered 
state of quantum antiferromagnets. Furthermore, we argue that vortices can 
be divided into two kinds: the first kind always carries zero angular momen- 
tum and are bosons, whereas the second kind carries angular momentum S 
under favourable conditions and are fermions if S is half- integer. A plausible 
consequence of our results relating to the RVB theories of High- Tc supercon- 
ductors is pointed out. 
PACS Nos, 75.10.Jm, 75.30.Ds, 79.60.-i 
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I. INTRODUCTION 



In the past few years there has been a lot of interests in the study of quantum anti- 
ferromagnets in two dimension on a square lattice, stimulated by the discovery of High-T^ 
superconductors Among others, one approximated approach to the quantum antiferro- 
magnet problems is the Schwinger-Boson Mean- field Theory (SBMFT). In this approach, 
the quantum spins are represented as Schwinger bosons and an approximated ground state 
is constructed by bose condensation of spin-singlet pairs. The mean-field theory can also 
be formulated in a large-N expansion as the saddle point solution to a generalized SU (N) 
quantum spin model [0,^. In the limit — oo, the mean-field theory becomes exact. 
SBMFT predicts that in one dimension, Heisenberg antiferromagnet is always disordered, 
with a non-zero spin gap in the excitation spectrum. The theory is found to give an adequate 
description for integer spin chains, while it fails to describe the massless state of half-integer 
spin chains. This is because the topological Berry phase term which plays crucial role in 
the latter case is not taken into account correctly in SBMFT. 0] However, in two dimension 
topological term does not exist and SBMFT is more reliable. In particular, the theory pre- 
dicts that the disordered (spin gap) phase exists only when the spin magnitude S is small 
enough, when S < Sc ^ 0.19. It turns out that SBMFT offers a fairly accurate description 
for the magnetic properties of the High-Tc compounds in the low doping regime ||^ . However, 
at large doping levels, the generalization of SBMFT which includes holes as fermions P|, 
does not seem to describe the High-Tc compounds correctly. For example, photo-emission 
experiments indicate that the compounds have "large" Fermi-surfaces which satisfy Lut- 
tinger theorem, whereas the generalization of SBMFT to include holes predicts a "small" 
Fermi surface with area proportional to concentration of holes. In the large doping regime, 
it turns out that the slave-boson mean-field theory |Q which treats the spins as fermions 
and holes as bosons produces a lot of properties of High-Tc compounds correctly |^ and 
the SBMFT is only superier in describing the low-doping, antiferromagnetic state. Thus a 
natural question is: what is the relation between the two theories? Can one understand the 
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two theories within a single framework? How does nature crossover from one description to 
the other as concentration of holes increases? 

In SBMFT, the spin operator Si on site i is expressed in terms of Schwinger bosons 
Si = ZiaZi, where a is the Pauli matrix, Zi{Zi) are two component spinors Zi = {Zi^, Zii), 
etc. Notice that in order to represent a spin with magnitude S, there should be 25' bosons 
per site 0]. The Hamiltonian can then be represented in terms of Schwinger bosons, and a 
mean-field theory can be formulated by introducing order parameters Ai^i^y =< ^«T^i+»7i ~ 
ZjjZj+^l > (?7 = ±x,±y) Alternatively, SBMFT can also be formulated in a large-N 
expansion as the saddle point solution to a generalized SU {N) quantum spin model. To 
derive the large-N theory, we divide the square lattice into A and B sublattices and consider 
the following transformation for the Schwinger bosons on i?-sublattice, 

for all sites j on the B sublattice. The Schwinger bosons on the A sublattice remains 
unchanged in the above transformation. The Lagrangian of the Heisenberg model can then 
be represented in the transformed boson coordinates as 

+ J E \^^,^+X-J E [^h+,Zf.Zf^,. + H.C.] (1) 

where Aj's are Lagrange multiplier fields enforcing the constraint that there are 25* bosons 
per site, Aj j+^'s are Hubbard- Stratonovich fields introduced in decoupling of H. a is the 
spin index. A SU (N) spin theory can be formulated with the above Lagrangian if the usual 
SU{2) spins t and [ are extented to the SU{N) case where N-spin species are introduced. 
SBMFT can be considered as the saddle point solution of the path integral over L which 
becomes exact in the limit — ^> oo ||^. For infinite system with no defects, the saddle 
point solution has position independent Aj j+.^'s and A's and the solution is formally very 
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similar to the BCS solution for superconductors, except that in present case the pairing 
objects are Schwinger bosons but not electrons. Another important difference is that in the 
present case, the pairing bosons are always located on different sublattices, whereas no such 
distinction is found in the case of superconductors. 

The similarity between SBMFT and BCS theory leads us to ask the natural question 
of whether vortex-like solutions can be found in SBMFT as in the case of superconductors. 
More precisely, one may ask the question of whether one can construct a stable solution in 
SBMFT where Aj_i_|_^ has a phase structure 

where 9{x) is a smooth function of x and has a singular point at ~ such that for 
distance |a;| >> 0, 9{x) —>■ 9{x) + 2tt if one moves vector x around a close loop C which 

— * 

enclosed the point x' — 0. As in the case of superconductors, the simplest way to address 
this question is to construct a continuum description for SBMFT (analogous to Ginsburg- 
Landau theory for superconductors) and study the possibility of having vortex solution 
in the continuum approximation. In the following we shall perform such a study in the 
disordered state in SBMFT where the Schwinger bosons are not bose-condensed and spin is 
a good quantum number. In realistic High-Tc cuprates long range antiferromagnetic order 
is destroyed by introduction of charge carriers (holes). We shall not consider complications 
introduced by holes here and shall assume simply a disordered magnetic state with realistic 
spin magnitudes which can be described by the spin gap phase of SBMFT. We shall examine 
the properties of vortices within the continuum description. In section II we shall derive the 
continuum theory. Based on the continuum equations, we shall first study the properties 
of a single, non-magnetic impurity in the disordered state of SBMFT, where we shall point 
out the important difference between perturbations which are symmetric to the A and B 
sublattices, and perturbations which distinguish the two sublattices. In section III, we shall 
study static vortex excitations in our model. Based on the continuum theory, we shall 
argue that vortices are stable topological excitations in the disordered state of 2D quantum 
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antiferromagnets. We shall then show that there exists two kind of vortices in the continuum 
theory, corresponding to vortices centered on the mid-point of a plaquette, and vortices 
centered on a lattice point. We shall show that the properties of these two kinds of vortices 
are very different, because of the different symmetry with respect to the two sublattices. The 
first kind of vortex which is symmetric to the two sublattices can carry only zero angular 
momentum and is a boson, whereas the second kind of vortex distinguishes between the two 
sublattices can carry angular momentum S under favourble conditions and is a fermion if S 
is half-integer. In section IV we shall concentrate ourselves at the case S = 1/2 which is the 
physical case of interests and shall examine properties of a "liquid" of fermionic vortices. 
Based on simple symmetry arguments, we shall argue that the effective theory of a liquid of 
fermionic vortices has precisely the same form as the slave-boson mean-field theory for spins 
in the undoped limit. The case with finite concentration of holes will also be discussed. The 
findings in this paper will be summarized in section V, where we shall discuss a plausible 
scenerio of how the system crossover from a state described by SBMFT to a state described 
by slave-boson mean field theory upon doping. 



To derive the continuum theory we following Read and Sachdev and consider the 
representation where our system has two sites per unit cell and introduce the 'uniform' and 
'staggered' components in the A and A fields, where 



II. CONTINUUM DESCRIPTION FOR SBMFT 



A, 



i[0(z±^)+g±,(z±^)]e 



,i[e(i±f)+A±^(i±|)] 



(2a) 



where q. 



-V - 



-V — 



Ay and in momentum space (for k ^ 0) 



(2b) 



A{k) = \[X^{k)-\^{k)]. 



(2c) 
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The field describes uniform antiferromagnetic correlations whereas q{x) field de- 
scribes spin dimerization (spin-Peierls) effects [0. 6{x) and are fields describing the 
corresponding uniform and staggered phase fluctuations, respectively. All fields are slowly 
varying on the scale of a unit cell. Notice that the reason why four independent (real) fields 
are needed to describe fluctuations in one unit cell is precisely because we have divided the 
system into A and B sublattices. In the case of usual superconductors such a distinction is 
not present and only two real fields (or one complex scalar field) enters into the Ginsburg- 
Landau equation. Notice that V^(a;) = d^O{x) and A^(a;) can be considered as U{1) gauge 
fields coupling to the Schwinger boson Z's. The 'uniform' gauge field V^{x) couples to the 
bosons on the two sublattices with same gauge charge, whereas the gauge charges for the 
'staggered' gauge field A^{x) are opposite on the two sublattices. A single vortex solution 
centered at x' = corresponds to a solution of SBMFT with boundary condition 



or in terms of the gauge field, there is a 'uniform' gauge flux of n/2 flux quantum 
passing through the origin, similar to the case of superconductors. Notice that in SBMFT 
where < Aj^j+^ >^ 0, the gauge symmetry associated with the 'uniform' gauge field V^{x) 
{Z"^ — > Z^e^^, Z^ Z^e^^, ^ — 2d^T) is broken, whereas the 'staggered' gauge 
symmetry (Z^ Z^e^ , Z^ Z^e'^^, A^, - d^V) remains intact in SBMFT. The 

existence of vortex solution is tied with the broken symmetry gauge field V^{x) as in usual 
superconductors. 

At distance much larger than lattice spacing, fluctuations associated with the 'uniform' 
variables and Qr are unimportant. Thus we shall neglect the V0(x) and V^r(a;) terms in 
the following and derive a continuum theory for the rest of the variables. Notice that the 
variable is kept in our continuum theory since the term is essential for studying of vortices. 
In the continuum limit, the Lagrangian L becomes 




- E \<P{x){l-\{V,{x)V^{x)))Z^{x)Z^{x)-\m^^ 



li=x,y 
1 

7B/^\ ryBl 



+ -q,[Z^{x)D,Z^{x) - Z^{x)D;Z^{x)] + H.C. 



+eAx)iz^ix)z^ix) + z^ix)z^ix)) 

+ AAx)iZ^ix)Z^ix) - Z^{x)Z^{x))} + \ Y^Mxf + q,{xf] - AS9,{x)^ , (3) 

where = + iA^^ and D* = — iA^. The continuum Lagrangian Eq. (|^) can be 
further simphfied by introducing new boson fields 

The TTo- field can be integrated out safely at large distance and low energy ||^, leaving 
effective Lagrangian for ifj^ field, 



Y^jdr j d'x {0(x) {\D;i;,\' + m{xy\i;^\') + iM^))-\D^ - q^D^)i^t{Dl + q,D;)i^„ 
+ + q.i^f] - 45 (20(x)(l - ^{V,{x)V^{x))) + mixf^ } (4) 



where 0(f)m(f)2 = 20^(f) -4[1 - (l/8)(V;,(f)\/^(f)]0(f) and and 0i(f) = (l/2){0^(f) + 
2[1 — {l/8){V^{x)V^{x)](f){x)}. In the limit q^ = and A^ = 0, and neglecting gradient 
terms V0, V6'^, the above Lagrangian can be easily diagonalized resulting in an effective 
Lagrangian Lejj in terms of the (f){x) and V^(x) fields. The dynamical effect of the remaining 
terms can be obtained by looking at Gaussian fluctuations of the flelds around the saddle 
point solution A^, q^ = 0. ^ In particular, in the small m(x) ~ m — > limit, we obtain (see 
Appendix) 

Lpff = 
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jdrj fx |(o - 4(2S + m + + ^)m' 



+^Plu + ^cF/^r^M + 0(m2)| , (5) 

where ~ m, a(< 4(25* + and c are constants of order 0(1). The precise values 
of a, 6, c depend on the underlying lattice structure and cannot be obtained in a continuum 
theory. We have consider the realistic case = 2 in deriving the above expression. Notice 
that similar effective Lagrangian has been obtained by Read and Sachdev previously in 
studying effect of instantons The only difference here is that the V^V^ term which was 
not considered by Read and Sachdev is now retained. Recall that is the 'uniform' f/(l) 
gauge field coupling to the bosons and the V^jV^ term in Eq. just represents the Meissner 
effect associated with nonzero order parameter A in SBMFT. Notice that as in case of usual 
superconductors, we have choosen the London gauge V .V = in our derivation. The term 
can also be written in a gauge invariant way by replacing by the gauge invariant object 

— 2d^T in Eq. (^. This term does not contribute to the instanton effect discussed by 
Read and Sachdev but is crucial to the study of vortices. 

To understand more bout the continuum theory we first consider the properties of a single 
non-magnetic impurity in the disordered state of SBMFT. We shall assume that a non- 



magnetic impurity simply replaces a spin at site i by a non-magnetic object. The simplest 
way to model this in our formalism is to is replace the constraint equation ZZ = 2S on 
site i by ZZ = 25(1 — rii) = 0, where rii is the number of non-magnetic impurity on site i. 
It is than easy to see that the effect of nonzero rii in our effective Lagrangian Eq. (^) is to 
introduce in the Coulomb gauge an extra term —{2Sniei)Ar, where Cj = ±1, depending on 
whether the impurity is located on the A- or B- sublattices, [|1^ i.e., non- magnetic impurities 



appear as effective (staggered) gauge charges of magnitude 256 localized at the impurity 
sites. A corresponding electrostatic potential V{r) ~ {2Sei)ln{r /C,) (where ^ ~ m~^) will be 
induced around the impurity. To lower the electrostatic energy, 25 bosons will be nucleated 
out from the vacuum to screen the electric field, resulting in formation of local magnetic 



moment of magnitude S around a non- magnetic impurity ||T0|. Notice that this effect has 



been observed experimentally in the High-Tc compounds in the underdoped region upon 
substitution of Cu by Zn in the conducting planes |TT| . 

This result can also be understood in an alternative way by observing that the same 
physical effect should have occurred if instead of replacing a spin on site i by an non-magnetic 
impurity, we let the Heisenberg coupling Jj to go to zero for those bonds joining to site i. 
In this case the 'impurity' site i behaves as an non-magnetic impurity as far as the rest of 
the system is concerned. The only difference is that a free spin of magnitude S now remains 
on site i. However, the effect of this approach on our effective Lagrangian Eq. @ looks 
rather different compared with our previous approach. Instead of introducing a 'staggered' 
electric charge of magnitude 2Se on site i, a boundary condition Aj = 0(?7 = ±x, ±y) is 
now imposed on our system. Using Eq. (pa]), the corresponding boundary condition in the 
continuum theory becomes 

where z — in the continuum limit, n is a unit vector. ± depends on whether site i is 
on A- or B- sublattices. The properties of the system around the non-magnetic impurity is 
obtained by minimizing the free energy Eq. under this boundary condition. 

Performing the calculation, we find first of all that 7^ and correspondingly q{xi) 7^ 

0. The reason for this behaviour can be understood easily from Eq. @ by noticing that 
the coefficient in front of the g^g'^ term becomes negative as becomes smaller than 2b, 
implying that spontaneously dimerization will occur when becomes small enough. Because 
of stability requirement, it can be shown that dimmerization |g| cannot have value larger 
than qo = (\/2)0 (see Appendix). Anyway both (p{xi) and q{xi) ~ should be nonzero 
in the continuum solution. The reason why q should be nonzero can also be understood 
from symmetry consideration. Notice that the (f){x) field is symmetric under exchange of A- 
and B- sublattices whereas g^(x) field is antisymmetric. Thus a solution with qfj,{x) field 
identically zero would not distinguish between the two sublattices. However, the single non- 
magnetic impurity problem certainly distinguishes the two sublattices since the impurity 
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can only be placed on either one of the two sublattices. Thus we expect to be nonzero 

around a single impurity. Minimizing the free energy with respect to the field we find also 
that E{x) ~ ^(x), where = iF^r is the 'staggered' electric field in the system. With the 
boundary condition q ~ (±r?,)0 around the impurity we find that electric field radiating out 
(or in) from the impurity is obtained in our solution, i.e., depending on whether the impurity 
is located on the A- or B- sublattices, it behaves as source or sink for the electric field, in 
exact agreement with our previous approach which predicts that non-magnetic impurity 
behaves as a single staggered gauge charge added to the system. 

One may also extend our approach to discuss the case when the Heisenberg coupling 
J' joining to the site i is small but nonzero, so that the spin at site i still couples weakly 
with the surrounding environment. The previous discussions still apply, except that the 
spin at site i appears now as a localized spinon state occupied by 25* spinous in the system 
(notice that at two dimension, an arbitrary small attractive interaction is enough to generate 
a bound state ^^). The spinon state again behaves as a (staggered) gauge charge (with 
magnitude 25"), and will bind 25* other spinous with opposite gauge charge (or localized on 
opposite sublattice) next to itself as in the case with J' = 0. However for nonzero J' the 
localized spinous on the two sublattices will interact antiferromagnetically because of the 
underlying antiferromagnetic correlation in the system, [|l^ forming at the end a localized 



spin singlet around site i. It is only in the limit J' ^ that the two spinous are decoupled 
and free local magnetic moments appear. 

Our discussion can also be extended to study other forms of local defects. For example, 
one may study the situation when the Heisenberg coupling Ji^i+x between sites i and i + x is 
set to zero, i.e., a single bond is being removed from the system. Notice that the perturbation 
is symmetric with respect to interchange of A- and B- sublattices. It is straightforward to 
generalize our previous discussion to this case where we find that the perturbation now 
appears as an effective electric dipole moment in Lejj. No local moment is expected to form 
around the impurity bond in this case because of the follow reason: because of the symmetry 
between A- and B- sublattices bosons(spinons) can be nucleated out from vaccuum only in 
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pairs, with one localized on the A sublattice and one on the B sublattice. The interaction 
between the nucleated spinon pairs will again be antiferromagnetic and has a magnitude 
~ Je~^^^, where x is the distance between the two spinous and ^ ~ m"^ is the 'size' of the 
spinon wavefunction. For the single bond defect, the distance x between the two spinous is 
of order of lattice spacing << ^. Thus the effective interaction between the two spinous is 
of order J, i.e., they will form a strong local spin singlet and no isolated magnetic moment 
will appear. 

Experimentally, the Cu ion in the conducting plane of high- Tc cupartes can be substituted 
by Zn (non-magnetic impurity) and magnitude of Heisenberg exchange J can be modified 
locally by introducing impurities out of conduction planes or substitution of O ions. It 
is observed that while substitution of Cu by Zn in the conducting plane introduces local 



magnetic moments in the underdoped (spin gap) phases of High-Tc cuprates |jTT[, no such 
effect is observed in other ways of introducing defects. Our theoretical investigation on 
disordered state of quantum antiferromagnets based on continuum description of SBMFT is 
in satisfying agreement with experimental results. ITU] 



III. STATIC VORTICES IN DISORDERED STATE OF SBMFT 

In the continuum theory, a unit vortex centered as a; = corresponds to a solution 
of SLgff = 0, with 6{r, Q) = Q, where r and Q are the distance and angle in the (2D) 
polar co-ordinate. The nonzero vorticity introduces a diverging kinetic energy through the 
V^V^" ~ (V6')2 term in Lejj, which supresses the magnitude of (j) around the vortex core 
r — i> 0, as in the case of usual superconductors. Minimizing the free energy with respect 
to 0, and keeping in mind that for small 0, |g| ~ (^2)0, we find that = |g| = for 
r < To ~ [(8(25* -|- 1) -|- 46 — 2a) /{2S + 1)]^^^, but not going to zero smoothly as r ^ 0, 
as in the case of usual superconductors. This result is an artifact of our continuum theory 
where the (V0)^ term is not included in Lejj. Nevertheless, the qualitative effect where 
is being suppressed around the vortex core is clear. For large r, 0(r) — ^ 0o = (4(25" + 1) — 
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a)/2 and g — > 0, indicating that vortices are stable topological excitations in SBMFT. The 
behaviour of vortices in SBMFT is very similar to vortices in usual superconductors at large 
r. However, the small r behaviours are very different. The suppression of (j) and q fields 
around vortex core reflects the fact that around the center of vortex, the bond amplitude 
Aj j+^'s are being suppressed. As has been discussed in the previous section, suppression 
of local bond amplitudes may lead to generation of localized magnetic moments, depending 
on the detailed bond configuration. To understand the properties of a vortex, it is thus 
important to understand the underlying bond structure at the vortex core. 

Before studying the bond structure at vortex core, let us explain first how local magnetic 
moment binding to vortex core modifies the properties of a vortex. For a unit vortex the 
magnetic flux seen by a Schwinger boson is half flux quantum, (recall that is an 'uniform' 
gauge field which has same gauge charge for bosons on both sublattices) Therefore the orbital 
angular momentum of a Schwinger boson around a unit vortex is ^ + integer. The composite 
object of one Schwinger boson bound to the vortex is thus a fermion because of the angular- 
momentum-statistics relation. Notice that this phenomenon has a direct correspondence 



in superconductors where the bound state of an electron to a vortex is a boson. However, 
if instead of one boson, two bosons are bound to a vortex, then the composite object is again 
a boson since angular momentum will be of integer value again. Now it should be clear why 
presence of local magnetic moment at vortex core is an importance issue. The magnitude of 
the magnetic moment measures the number of Schwinger bosons bound to the vortex and 
determines the statistics of the vortex. 

The bond structure at the vortex core depends on where the center of vortex is being 
located. Here we shall consider only vortices which have maximum rotation symmetry 
allowed by the underlying lattice, i.e. rotation by n7r/2. In this case there are only two 
possible type of vortices, with vortex center located (i) at center of a plaquette and (ii) at 
a lattice point, (see Fig.l) We shall first consider vortex (i). 

From continuum theory it is expected that the four bonds surrounding the center of 
vortex (see Fig. la) will be largely suppressed because of the kinetic energy associated with 
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the V/jV'^ term. Following similar analysis as in previous section we find that the sup- 
pressed bonds introduce in the continuum theory an effective electric quadrapole structure 
surrounding the vortex center. Notice that as in the case of single broken bond the structure 
is symmetric with respect to the A- and B- sublattices. Thus we expect that localized spinon 
pairs may form around the center of vortex but isolated magnetic moments cannot occur. 
Thus this kind of vortex carries zero angular momentum and is a boson. 

The second kind of vortex is centered at a lattice point and as a result, it is expected 
that the four bonds joining to the vortex center (see Fig. lb) will be largely suppressed. 
This situation is very similar to the case of single non-magnetic impurity discussed in the 
last section. First of all, 25" spinous will be found localized at the vortex center because of 
the suppressed bonds. The spinous behaves as a 'staggered' gauge charge with magnitide 
25" and 2S other spinous with opposite gauge charge will be nucleated from the vacuum 
to screen the 'staggered' electric field generated by the spinous localized at the vortex cen- 
ter. The spinous on opposite sublattices interact antiferromagnetically through an effective 
Heisenberg exchange of order J' ~ J| Acp, where A^, is the bond amplitude of the four bonds 
joining to the vortex center. As a result, a spin singlet will be formed. The vortex carries 
zero angular momentum and is again a boson. Notice that this is true even in the limit 
J' — > since integer number of bosons = 43 are found binding to the vortex. 

The situation is however quite different if we consider a finite density of the second kind 
of vortices. We shall show that when density of vortices is large enough, it may become 
energetically unfavourable to nucleate spinous from vacuum to screen the staggered gauge 
field and as a result, only 25" spinous will be found binding to vortex center resulting in 
fermionic vortices when S is half-integer (odd number of bound bosons) . For simplicity we 
shall first consider two vortices seperated by distance /, with one vortex on each sublattice. 
First let us consider the energy when spinous are nucleated from the vacuum to screen the 
staggered gauge charge on each vortex. The total energy will be sum of three terms: (i) 
the electrostatic energy which is of order 2(2Se)'^ln{lo/^), where ~ ^ is the 'size' of the 
nucleated spinon wavefunction, (ii) The exchange energy between the spinous localized as 
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vortex center and nucleated spinons, which is of order —2S{S + 1) J'e~'°/^ ~ —2S{S + 1) J', 
and (iii) The energy needed to nucleate the spinons from the vacuum, which is of order 
iSm. The sum of the three terms is of order 2S{2m — {S + 1) J')- The energy for the second 
case when no spinons are nucleated from vacuum consists of two terms: (i) the electrostatics 
energy which is of order {2Se)^ln{l/$,) and (ii)the exchange energy between spinons located 
at the center of the two vortices which is or order —S{S+ 1) J'e~'/^. In 2D, ~ m and the 
electrostatic energy is of order AS'^m.lnQ,/^). For I » ^, it is certainly energtically more 
favourable to nucleate bosons from vacuum to screen the vortex gauge charge. However, 
when / < 4, the energy cost of the second case is of order —S{S + 1) J', and is energetically 
more favourable if S{S + 1)J' < ASm. Notice that J' ~ -^^l^cP is expected to be very 
small because of suppresion effect around vortex core. Thus the condition can be satisfied 
easily even with a relatively small m. For finite density of vortices S, I ~ 6^ Thus for 
low concentration of vortices, we expect that spinons will be nucleated from vacuum to 
screen the staggered gauge field and the second kind of vortices are bosons. However, when 
density is high enough {5 ~ (1/^)^) and m large enough, the situation will change and it 
may be energetically more favourable for the vortices to stay as fermions (for half-integer 
S). Notice that finite density analysis for the first kind of vortices does not show similar 
qualitative difference between low and high density. Because of the symmetry between A- 
and B- sublattice, bosons are always added or removed from a single vortex in pairs which 
cannot affect the statistics of a vortex. 

Before ending this section let's summarize our findings here. Within a continuum de- 
scription we have demonstrated in SBMFT the existence of vortex solutions which are stable 
topological excitations. We have found two different kind of vortices, according to the dif- 
ferent locations of vortex center. The first kind of vortex has center located at center of a 
plaquette whereas the second kind of vortex has center located at a lattice point. The two 
kind of vortices have very different characters. The first kind of vortex is always a boson, 
whereas the second kind of vortex can be a fermion (for half-integer S) under conditions 
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which are rather easy to acheive. We cannot, however, be absolutely sure whether 'fermionic' 
vortices exist because of the natural limitation of a continuum theory which gives only or- 
der of magnitude estimates. Notice also one important distinction between the vortices we 
study in this paper and usual vortices in superfluids. In the present case, vortices cannot 
be distingished by sign of vorticity since vortices carrying 1/2 (or tt) and —1/2 flux quanta 
should be considered as identical. On a lattice where the order parameter field Aj is a 
link-variable, vr and —tt vortices can be related by a pure gauge transformation. The usual 
singularity encountered in the gauge transformation in continous space does not arise here 
because of the discretized lattice structure. 

It is important to clarify the different roles played by the 'uniform' and 'staggered' 
gauge fields in deciding the properties of vortices. The gauge fields arise from fluctuation in 
phases of the order parameter Aj j+^'s. The uniform gauge field couples to bosons on the 
two sublattices with same gauge charge, and the corresponding gauge symmetry is broken 
in SBMFT. Vortices are stable because of this broken gauge symmetry, and corresponds 
to solutions with "uniform" magnetic flux penetrating center of vortices. The 'staggered' 
gauge field couples to bosons on opposite sublattices with opposite gauge charge. It plays 
no role in generating a stable vortex solution, but has strong effects in determining the 
precise properties of the vortex core. Without the 'staggered' gauge field, spinous need not 
be nucleated from vacuum to screen the 'staggered' charges generated at, for example, the 
center of the second kind of vortex. In this case, for half-integer spin systems, the second 
kind of vortices will always be fermions independent of density, since only odd number of 
bosons are found binding to center of each vortex (corresponding to the spin S at center of 
vortex). Similarly, local magnetic moments will not be generated by non-magnetic impurities 
if 'staggered' gauge field does not exist. 

It is also interesting to point out that similar vortex excitations have been considered 
by Read and Chakraborty ||T5[ in a short-ranged RVB wavefunction for 5 = 1/2 quantum 



antiferromagnets. They considered also the two kinds of vortices we discussed here. The 
statistics of the vortices were examined by direct Berry phase computations where similar 
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conclusion that the first kind of vortex is a boson, whereas the second kind of vortex is 
fermion independent of density is arrived at their study. The effect of 'staggered' gauge field 
on vortices is not included in their study. The vortices we discuss in the present paper can 
be viewed as generalization of their vortices in short-ranged RVB wavefunction to the more 
complicated situation described by SBMFT. 



IV. QUANTUM LIQUID OF FERMIONIC VORTICES 

In this section we shall examine the properties of a quantum liquid of fermionic vortices 
based on general symmetry considerations. We shall not ask the dynamical question of how 
the liquid of vortices come to exist in the first place but shall ask the simpler question if the 
liquid exists, what general properties should it has? For simplicity we shall restrict ourselves 
to the case of S* = 1/2 antiferromagnets which is the physical case of interests. In this case, 
a unit fermionic vortex carries spin ±1/2. Notice that vorticity does not contribute another 
independent quantum number as discussed in last section. 

First of all we shall consider the Fock space for vortices. The 'vacuum' state for vortices 
is just a spin-disordered state in SBMFT. A one vortex state with spin a at position i is 
a vortex centered at site i with spin pointing in direction a. Notice that two vortices with 
same a seperated by distance d are not orthogonal to each other but have finite overlap 
< Xi\x2 >~ e~'^l^. To construct orthogonal set of vortex states, we consider Wannier vortex 
states which can be constructed from vortex states using standard methods. The Wannier 
states are orthonormal to each other and represent localized vortices given that the overlap 



< X\\x2 > falls off exponentially as function of distance, [pig 

Next we shall introduce second quantization representation for vortices. The creation 
operator cf^ is defined as an operator which turns on an 'uniform' magnetic flux of 1/2 flux 
quantum centered on site i, and creates a Schwinger boson with spin a on site i. Similarly, a 
destruction operator Qo- turns on an 'uniform' magnetic flux of 1/2 flux quantum in opposite 
direction, and removes a Schwinger boson with spin a on site i. With the constraint that we 
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cannot put more than one Schwinger boson on any site, it can be checked easily that these 
operators obey usual fermion commutation rules. Notice that because of the constraint that 
there are always one Schwinger boson per site, a physical vortex cannot be created by cf^ 
operating on the Vacuum' state, which changes the number of Schwinger bosons on site i. 
An example of an operator which creates a 'physical' vortex state from vacuum is cf^Zi^, 
where Zicr is a Schwinger boson destruction operator with spin a on site i. Notice also that 
the vortex occupation number < n^^ >=< cf^Cia > is equal to the corresponding Schwinger 
boson occupation number < n^^ >=< Zi^Zifj >, since the only difference between vortex 
and boson operators is the magnetic flux which is being cancelled in the occupation number 
operator. 

To proceed further we argue that the dynamics governing the quantum liquid of vortices 
must possess a local SU (2) symmetry, where the Lagrangian for vortices must be invariant 
under the transformation 



with + = 1. Notice that the same SU{2) symmetry was obtained in the Heisenberg 
model for 5=1/2 spins when the spins are expressed in terms of Fermion operators. |[T6|Jl7[| 
The reason why we expect the local SU{2) symmetry to exist can be seen most easily by 
considering a physical vortex state with spin a on site i and ask the question: what is the 
corresponding 'hole' or 'antiparticle' state for this vortex state? To construct the 'hole' state 
we operate Qo- on the 'vacuum' state, which creates a magnetic flux of 1/2 flux quantum 
centered at site i, and remove a spin a Schwinger boson, leaving site i empty. However, this 
is not a physically allowable state since the physically allowable state must has one boson 
per site. To construct a physical 'hole' state we must put back a Schwinger boson on site i. 
Now we may ask, which Schwinger boson shall we put back? It cannot be a spin a Schwinger 
boson, since putting back a spin a Schwinger boson just gives the original 'particle' vortex 
state. The only admissible choice is thus to put on site i a spin —a boson, which creates a 




aicf^ + (3iCii 



(6) 
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spin —a 'particle' vortex state. Thus we expect that the 'hole' state for a spin a vortex is a 
spin —a vortex, which gives rise to the local SU (2) symmetry stated above. Another related 
argument is to observe that our underlying system is a pure spin system. Thus the effective 
vortex Hamiltonian must be equivalent to some effective pure spin Hamiltonian. The spin 
operator when expressed in terms of fermions possesses the SU{2) symmetry stated above 
p!6| , [T7|| , so must be a pure spin Hamiltonian. 



The SU (2) local gauge invariance imposes a very severe constraint on the plausible form 
of the effective theory for the quantum liquid of vortices. A mean-field theory of a quantum 
liquid which is a spin singlet must describe a gas of fermions with particle-hole symmetry, 
i.e., a half- filled band with a particle- hole symmetric band-structure. Fluctuations around 
the mean-field theory must be described by a SU (2) gauge theory. Physical excitations are 
represented by particle-hole excitations (generation of vortex pairs). These properties are all 
basic ingredients of the slave-boson-mean-field theory at half-filling, when holes are absent. 
I6| , p!7[] Notice in particular that the mean-field theory predicts that the Fermi surface for 



the vortices, if exists, must be 'large', i.e., it obeys the Luttinger theorem at half- filling, 
since < n^^. >=< >= 1/2 for spin-singlet wavef unctions. 

In the presence of holes the discussion is similar if one assumes that holes form bound 
state with vortices. Notice that in SBMFT where spins are bosons, holes are fermions and 
can be represented by 'slave-fermion' operators and /j. A hole binding to a spin a 
vortex is formed by removing a spin a electron from the center of the spin a vortex. The 
resulting object is a boson, since it is formed by binding a fermion (hole) to a magnetic 
flux of 1/2 flux quantum. As in the case of vortices, we may introduce boson creation and 
destruction operators bf and bi for the vortex-hole bound states. A boson creation operator 
bf inserts a magnetic flux centered at site i, and create a (fermionic) hole on site i. Notice 
that as in the case of vortex operators, a physical hole cannot be created by operating bf on 
the vacuum alone without first removing the spin on site i. An example of an operator which 
creates a physical vortex- hole bound state from vacuum is thus bfZi^j, where the spin on 
site i is removed by the Z operator. Notice also that the vortex-hole bound state occupation 



number < bfbi > is equal to the corresponding fermionic hole occupation number < f^fi >, 
since the only difference between the vortex-hole bound state operator and hole operator 
is the appearance of the magnetic flux, which is again cancelled in the occupation number 
operator. 

In the presence of holes the SU{2) symmetry associated with vortices will be broken. 
To see that we repeat our construction of the 'hole' state for a spin a vortex. In presence 
of holes the state generated by Qo- operating on the vacuum is a physical state, since now 
holes are allowed in the system. Thus the 'hole' state for a spin a vortex is not necessarily 
a spin —a vortex, and the SU{2) symmetry is broken. The only symmetry which remains 
is the f/(l) symmetry Ci„ Ci„e^^,bf bfe^^^, which follows from the requirement that 
a physical electron destruction operator is the product of the spin-annihilation operator 
and hole creation operator and that all physical processes must be expressed in terms of 
electron operators. Notice that because of the no double occupancy constraint < n^^ > + < 
n^l > + < bfbi >= 1; a mean field theory which treats the vortices and vortex-hole bound 
states as independent gases of fermions and bosons, respectively will produce a vortex fermi 
surface (if exists) obeying the Luttinger theorem < n^^ >=< n^^ >= (1 — S)/2, where S 
is the concentration of holes in the system, as in the slave-boson mean field theory for the 
t - J model §. 

V. CONCLUSION 

In this paper we study the properties of vortices in disordered state of SBMFT. Based 
on the continuum equations we show that vortices are stable topological excitations in the 
system. Furthermore, we show that because of the underlying lattice structure, there exists 
two kind of vortices with vortex center located at (i) center of a plaquette, and (ii) lattice 
point. The two kind of vortices have different properties because of the different symmetry 
with respect to the two sublattices. In particular, for half-integer spin systems, we argue 
that the second kind of vortices can be fermions, under conditions which are not difficult 
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to achieve. We next consider a quantum liquid of fermionic vortices and show that because 
of general symmetry requirement, the effective theory for the liquid of vortices must has 
a general form very similar to the slave-boson mean field theory in the undoped limit. In 
the presence of holes similar conclusion is reached if one assumes that holes form bound 
states with vortices. A very interesting possibility brought out by our analysis is that the 
Schwinger-boson-slave-fermion mean field theory which is believed to describe the t — J 
model qualitatively correctly in the 5 — » limit and the slave-boson mean field theory which 
is believed to describe the metallic state more correctly can be understood within a single 
framework, if one identify the fermionic spins in the slave-boson mean field theory as vortices 
in Schwinger-boson treatment. In this paper we have demonstrate the possibility of such an 
identification. However, we have not been able to show within the framework of Schwinger- 
boson-slave-fermion mean field theory that fermionic vortices and vortex-hole bound states 
are more important low energy excitations then bosonic spinwaves and fermionic holes when 
density of holes is high enough. The problem is related to the more fundamental question of 
how the antiferromagnetic spin background get disordered when holes are being introduced. 
We shall not try to address these questions further in this paper, except to point out lastly 
that based on the slave-fermion mean field theory for the t — J model, fermionic vortices 
and vortex-hole bound states indeed become more stable in the presence of holes. 

In the presence of holes, the effective Hamiltonian for the spins (in Schwinger-boson- 
slave-fermion treatment) has a form 



where rj = ±x,±y,a = A,B,Q = < Zf^Z^_^ > and F =< ft fi+r, >• The first 



term is the original Heisenberg interaction and the tF + JQ term is an effective hopping 
term for spins generated by motion of holes. The second term represents the constraint 



H 



+AE E ^r.^r. - (1 - ^) 



a 




(7) 
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ZiZi + ff ji = 1 satisfied on average. Notice that the effective Hamiltonian for Schwinger 
bosons still has a form very similar to the BCS Hamiltonian for superconductors. Thus 
we expect stable vortex solutions will exist. The most important difference in the presence 
of holes is the existence of the effective hopping term which breaks the staggered gauge 
symmetry. [|T^ The 'staggered' gauge charge is no longer a conserved quantity and as a result, 
the 'staggered' gauge field we discuss in the previous sections will no longer be important 
TSj . In particular, the second kind of vortex will be more likely to stay as fermions since 



the 'staggered' electric field which is responsible for binding another boson to the vortex is 
no longer important. It is also easy to see that provided that fermionic vortices exist, it is 
energetically favourable to bind holes to them. The reason is simple: the antiferromagnetic 
correlation around the vortex core is suppressed by the V^V^ term and as a result, a hole 
can gain more kinetic energy by staying at the center of a vortex. 

This work is supported by Hong Kong UGC Grant HKUST636/94P. 



APPENDIX A: 

In this appendix we outline our derivation leading to the effective Lagrangian Eq. (^. 
First we shall neglect the staggered gauge field and shall set the terms q^{x) and 

V^(x) equal to constants, i.e. we neglect gradient effects associated with these terms. In 
this case the effective Lagrangian Eq. @ can be diagonalized easily leading to Free energy 
(at zero temperature) 

F = 2Y,E(k) + ^(202 + g.g) - 2(25 + 1)(20(1 - \v.V)+m'') (Al) 
^ z o 

k 

where E{k) = +{2(f)'^k.k — {q.ky+m'^y^'^ is the energy for the (bosonic) spinous, m is the spin 
gap. The factor 2 is coming from spin sum. In the limit m and g — > 0, E{k) y/2(f)\k\ 
and F {a- 4(23 + 1))0 + 0^ + (5 + l/2)0(F^y), where a is a cutoff dependent constant 
coming from J2E{k). For q small, we can expand the {q.kY term in E{k) to obtain the 
q^^q^ term in Eq. (^. Notice that is always a postive real number since the postive root 
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— # — * 

is choosen in computing E{k). Notice also that the stabihty criteria E{k) > for all values 
of k leads to the condition q.q < 20^. The terms involving the staggered gauge field in 
Eq. (|^) can be obtained by considering Gaussian fluctuations associated with the field. 
Detailed calculation can be found in reference |^ and we shall not repeat the derivation 
here. 
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FIGURES 

FIG. 1. two possible kind of vortices with vortex center located at (a) center of a plaquette, 
and (b)latticc point. In (a), the four bonds surrounding the vortex center and in (b), the four 
bonds joining to the vortex center, are suppressed. 
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